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Abstract. We present the results of our recent analysis of the meson-photon transition form factors 
Fpy(Q 2 ) for the pseudoscalar mesons P = n°, rj , rj' , rj c , using the local-duality version of QCD sum 
rules. 
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^ INTRODUCTION 

The processes 7* 7* — > P with P = 7T°, r/ , r/ ; , r\ c are of great interest for our understanding 
of QCD and of the meson structure. In recent years, extensive experimental information 
on these processes has become available [1-5]. 

The corresponding amplitude contains only one form factor, FpyJq\,q^)\ 

. 

g: (f(qi)f(q2)\P(p)}=i£e,e 2qiq2 F PYr (qW 2 ). (1) 

1 1 

A QCD factorization theorem predicts this form factor at asymptotically large spacelike 
momentum transfers q\ = — Q\ < 0, q\ = — Q\ < [6]: 

F Pn (QlQ 2 2 )^2el J Q 2^2 ( ( f^ r # SS (^) = 6K(1-^). (2) 

X: 1 2 

Hereafter, we use the notation Q 2 = Q 2 , and < /3 =Q\/Q\<1 (that is, Q 2 is the larger 
virtuality). For the experimentally relevant kinematics Q 2 rj and Q 2 = Q 2 , for instance, 
the pion-photon transition form factor takes the form 

Q 2 Fny(Q 2 ) -> V^A, A = 0.130 GeV. (3) 
Similar relations arise for r\ and r\' after taking into account the effects of meson mixing. 

DISPERSIVE SUM RULES FOR THE ff^P FORM FACTOR 

The starting point for a QCD sum-rule analysis of the 7* 7* — > P transition form factor is 
the amplitude 

(0\jl\f(q2)f(qi)) = e% aj5 (p\q h q 2 )£?£ 2 3 , P = qi+qi, (4) 



where £12 are the relevant photon polarization vectors. This amplitude is considered for 
— q\ = Q 2 > and — q 2 = Q\>0. Its general decomposition contains four independent 
Lorentz structures (see e.g. Refs. [7, 8]) but for our purpose only one structure is needed: 

Tfta^ipWuqi) =PLi £ ap qiq2 iF(p 2 ,Qi (£) + ■■■. (5) 

The corresponding invariant amplitude F(p 2 , Q 2 V Q 2 ,) satisfies the spectral representation 
in p 2 at fixed Q 2 and Q 2 

CO 

F{p\QlQ 2 2 ) = - [ -^A(s,Q 2 ,Q 2 ), (6) 

•?th 

where A(s, Q\, Q 2 ) is the physical spectral density and s± denotes the physical threshold. 
Perturbation theory yields the spectral density as a series expansion in powers of a s : 

A p qcd(*, <2i, Ql\m) = A^ CD (5, Qj, Q\\m) + *A$a>(s, Ql Q 2 2 \m) + ■ ■ ■ , (7) 

where m is the mass of the quark propagating in the loop. The lowest-order contribution, 
Apqcd( 5 ' Qv Q\\ m )' corresponding to a one-loop triangle diagram with one axial current 
and two vector currents at the vertices, is well-known [9]. The two-loop 0{a s ) correction 
to the spectral density was found to vanish [10]. Higher-order corrections are unknown. 

The physical spectral density differs dramatically from A p qcd(s, Q 2 , Q 2 ) in the low-s 
region; it contains the meson pole and the hadronic continuum. For instance, in the / = 1 
channel, one has 

A(s, Ql Ql) = k8(s- ml) V2f K F nn (Ql Q 2 ) + 9(s- s th ) a£J(*, Ql Ql) ■ (8) 

The method of QCD sum rules allows one to relate the properties of the ground states to 
the spectral densities of QCD correlators. The following steps are conventional within 
the QCD sum-rule method [11, 12]: equate the QCD and the physical representations for 
F(p 2 , Q 2 , Q 2 ,); then perform the Borel transform p 2 — > T, which suppresses the hadronic 
continuum; in order to kill then potentially dangerous nonperturbative power corrections 
which may rise with Q 2 , take the local-duality (LD) limit x = [13]; finally, implement 
quark-hadron duality in a standard way as low-energy cut on the spectral representation, 
in order to arrive at the following expression for the ground-state transition form factor: 

Kf P F Pyj {QlQ 2 2 )= J dsA pQCD (s,QlQj\m). (9) 

4m 2 

All details of the nonperturbative-QCD dynamics are contained in the effective threshold 
s e ff(<2j, Q 2 ). The formulation of reliable criteria for fixing effective thresholds proves to 
be highly nontrivial [11]. 

At large Q\ = Q 2 ->• <*> and fixed ratio ^ = Q\ /Q 2 , the effective threshold s eff (Qj , Q\) 
may be determined by suitable matching to the asymptotic pQCD factorization formula. 



From this, one finds that, in the general case m ^ 0, s s ff{Q 2 — > °°, j8) depends on /3. The 
only exception to this is the case of massless fermions, m = 0: in this case the asymptotic 
factorization formula is reproduced for any /3 if one sets 5 e ff(2 2 = 4tt 2 /|. The 

LD model for the transition form factor emerges when one assumes that, at finite values 
of Q 2 , s e ff(Q 2 ,j5) may be sufficiently well approximated by its value for Q 2 — > <», that is, 

Seff(£ 2 ,/3)=*eff(0 2 ^~,/3). (10) 

Introducing the abbreviation Fpy(Q 2 ) = Fpyy{0, Q 2 ) for the pseudoscalar-meson-photon 
transition form factor, its LD expression for Q 2 = and m — reads, in the single-flavour 
case, 

Independently of the behaviour of s e ff (Q 2 ) at Q 2 — > 0, Fp Y (Q 2 = 0) is related to the axial 
anomaly [7]. 



THE TRANSITION yf P IN QUANTUM MECHANICS 

The accuracy of the LD model for the effective threshold may be estimated in quantum 
mechanics. There, the form factor may be found exactly by some numerical solution [14] 
of the Schrodinger equation. From this, the exact effective threshold may be calculated: 
for any given experimental or theoretical form factor, the corresponding exact effective 
threshold is defined as the quantity that reproduces this form factor by a LD sum rule (9). 

The result from a quantum-mechanical model with a harmonic-oscillator potential [7] 
is shown in Fig. 1. For "light" quarks, the LD threshold gives a very good approximation 
to the exact threshold for Q > 1-1.5 GeV. For "charm" quarks, the local-duality model 
works for Q > 2-3 GeV. The accuracy of the LD approximation further increases with Q 
in this region. 
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FIGURE 1. The exact effective threshold in quantum mechanics, k^{Q), for two different values of the 
nonrelativistic constituent quark mass rag. 



ff^T] c FORM FACTOR 



In the case of massive quarks, we may exploit not only the correlation function (AW) as 
in Eq. (4) but also the correlation function (PVV) [8]. For each of these objects, an LD 
model may be constructed. By matching to the pQCD factorization formula, we derive 
Seff(Q 2 — > °°, j3) f° r (AW) and for (PVV). The results of the corresponding calculation 
for rj c are depicted in Fig. 2. Obviously, the exact effective thresholds corresponding to 
(AW) and (PVV), s^ v (Q 2 — >■ °°,/3) and s^ v (Q 2 °°,/3), differ from each other; they 
also differ from the effective thresholds of the relevant two-point correlation functions. 

Assuming that 5 e ff(2 2 ,/3) = Seff(<2 2 — °°,/3), we obtain the results shown in Fig. 2. For 
the above reasons, at very small Q 2 the applicability of our LD model is not guaranteed. 
Nevertheless, applying our LD model down to Q 2 = predicts F-q c y(Qi) = 0.067 GeV 1 
from the analysis of (AW) and P^O) = 0.086 GeV 1 from the analysis of (PVV); this 
has to be compared with the experimental number Pr; c y(<2 2 = 0) = 0.08 ± 0.01 GeV 1 . 
Seemingly, the LD model based on the correlator (PVV) gives reliable predictions for a 
broad range of momentum transfers Q 2 starting even at very low values of Q 2 (cf. [15]). 
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FIGURE 2. Form factor for the transition 77* — > r] c : exact effective thresholds s^ v (Q 2 — > °°, j3 ) (a) and 
s es V (Q 2 °°j P ) (b); form factors obtained for finite Q 2 from the LD sum rules for the correlators (AW) 
and (PVV) (c); LD model for the correlator (PVV) confronted with experimental data by BaBar [3] (d). 



77* 07 ^0 FORM FACTORS 

Here, the mixing of strange and nonstrange components [16] must be taken into account: 
Ffjy = F n y cos (j) —F s y sin (j), Fj.iy = F n y sin (j) +F s y cos (j), (j) ~ 38°, (12) 
with n — > (uu + dd) / a/2 and s — >• ss. The LD expressions for these two form factors read 



F„ y (!2 2 ) = i / dsA„(s,<2 2 ), F sy (Q 2 ) = i / dsA s (s,Q 2 ). (13) 



Accordingly, two separate effective thresholds emerge: s^j = 4n 2 f 2 , = 4n 2 f 2 , with 
/„ ps 1.07/n, / 5 ~ 1.36/^. The outcomes from the LD model [7, 8] and the experimental 
data [1, 4] are in reasonable agreement with each other (Fig. 3). 
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FIGURE 3. 77* — > (TJ,Tj') transition form factors P{n,r\')y{Q^)'- LD predictions [7, 8] (dashed lines) and 
recent fits [17] (solid lines) to the experimental data [1,4]. 



77* 71° FORM FACTOR 

First of all, we emphasize that the large- Q 2 behaviour of the r\, J]', and 7t° form factors is 
determined by the spectral densities of perturbative QCD diagrams and should therefore 
be the same for all light pseudoscalars [17]. In order to demonstrate this, we observe that 
the sum rule for (AW) in the LD limit % = is equivalent to the anomaly sum rule [18] 



Fny(Q 2 



1 



2V2n 2 f„ 



1 -2n 



Q 



(14) 



Similar relations arise for the 1 = and the ss channels. As shown in Ref. [17], the form 
factors F K y(Q 2 ), F r] y{Q 2 ), and iy y (<2 2 ) at large Q 2 are determined by the behaviour of 
the appropriate A con t(s, Q 2 ) at large s. By quark-hadron duality, the latter are equal to the 
corresponding A p qcd(s, <2 2 ); these are purely perturbative quantities and therefore equal 
to each other for the different channels. 
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FIGURE 4. 77* -> 7T° transition form factor F K7 (Q 2 ): LD prediction (magenta lines) and a fit [17] (solid 
lines) to the data [1, 2, 5]. The equivalent effective threshold s e f((Q 2 ) for each data point is found via (11). 



However, the B AB AR data for the % transition form factor exhibit a clear disagreement 
with both the r\,r\' form factors and the LD model at Q 2 as large as 40 GeV 2 . Moreover, 
opposite to findings in quantum mechanics, the violations of LD rise with Q 2 even in the 
region Q 2 m 40 GeV 2 ! We thus conclude that the B AB AR results are hard to understand 
in QCD (see also [19]). Noteworthy, recent Belle measurements of the 7ry form factor — 
although being statistically consistent with the BaBar findings (see [20, 21]) — are fully 
compatible with the r\ and r\' data as well as with the onset of the LD regime already in 
the region Q 2 > 5-10 GeV 2 , in full agreement with our quantum-mechanical experience. 



CONCLUSIONS 

We studied the 7t , r], r}', and r) c transition form factors by QCD sum rules in LD limit; 
the key parameter — the effective continuum threshold — was determined by matching the 
LD form factors to QCD factorization formulas. Our main conclusions are the following: 

• For all P — > yy* form factors studied, the LD model should work well in a region of Q 2 
larger than a few GeV 2 : the LD model works reasonably well for the 77 — > yy* , r\' — » yy*, 
and r\ c — > yy* form factors. For n° — > yy*, the BaBar data indicate an extreme violation 
of local duality, prompting a linearly rising (instead of a constant) effective threshold. In 
contrast to this, the Belle data exhibit an agreement with the predictions of the LD model. 

• Nevertheless, a better fit to the full set of the meson-photon form-factor data seems to 
prefer a small logarithmic rise of Q 2 F P j(Q 2 ) [17]. If established experimentally, this rise 



would require the presence of a l/s duality-violating term in the ratio of the hadron and 
the QCD spectral densities. 

• A high accuracy of the LD model has implications for the pion's elastic form factor: we 
can show that the accuracy of the LD model for the elastic form factor increases with Q 
in the region Q 2 m 4-8 GeV 2 [7]. The accurate data on the pion form factor suggest that 
the LD limit for the effective threshold, s e ff (<2 2 — > 00 ) = 4tt 2 /^, may be reached already 
at Q 2 = 5-6 GeV 2 . This property should be testable with the JLab upgrade CLAS12. 
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